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Abstract 


High voltage power transmission lines are important 
elements in an electric power system. Unavoidable accidents 
often occur resulting in faults on these transmission lines. 
During the past few years, a considerable amount of studies 
and many research projects have been undertaken to develop a 


digital transmission line fault location algorithm. 


A method involving curve fitting techniques for 
locating transmission line faults has been studied and 
tested in this project. The method is based on the least 
error squares approach which utilizes information of the 
first few cycles of post-fault data. The effects on the 
accuracy of the proposed fault location, due to the presence 
of line shunt reactance and fault resistance, have been 
investigated. Also, different order digital filters are 
employed in the testing to study the delay associated with 


each filter. 


A Gigital computer program has been successfully 
developed and tested on locating faults of digitally 
Simulated transmission lines. A fault detection routine is 
Incorporateqdse into. the; sprogram= to© initialize- the fault 


locating routine. 
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Cc - capacitance of transmission line in Farad 
Vv - instantaneous voltage 

i - instantaneous current 

V - voltage vector or phasor 

I Se CUGLeENCERVeCtOL, Or, phasor 

t - time in second 

y - length in mile 
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Superscripts 


Subscripts 


eS, 8 | 


pseudoinverse of matrix A 
inverse of matrix A 

error symbol 

summation symbol 

partial differentiation symbol 
time constant 


fundamental angular frequency in 


radian/seconda 


phase angle in radian 
factorial symbol 


Square root 


transpose of matrix 
Current quantity 
voltage quantity 
pre-fault quantity 


post-fault quantity 


kth component 

peak quantity 

line quantity 

fault quantity 

Quantity corresponds to total number 


line parameter corresponds to length l 
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- phase a to neutral quantity 
- phase b to neutral quantity 


- phase c to neutral quantity 


CHAPTER I 


Introduction 


A modern electric power system is” Fa ycomplex 
interconnection of generators, transformers, transmission 
lines, static and rotating loads, protective relays, circuit 
breakers and various control elements. Power can be 
transmitted either by direct current, single-phase or 
poly-phase alternating current. In practice, three-phase 
Sinusoidal alternating current has been adopted as the 
method of transmission due to its efficiency in generating 


power and in its transfer capacities. 


Electrical energy is a fundamental essence to the human 


society. A continuous supply of electricity, irrespective 
ClmeivitonmentalsCcONdItions, sis. amportantesito our. daily 
lives. During the occurrence of faults, power utilities 


must be able to isolate the faulty section from the rest of 
the system to ensure the continuity of electric power 
Supply. At the same time, a large current could result and 
flow through the network, which would cause severe damage to 
the system components. Damage to equipment and disruption 
of power could be minimised by the rapid disconnection of 
the faulty section by means of circuit breakers and 
protective relays. The type of relay or relays used in the 
protective scheme 1s mainly determined bythe nature of the 


protected equipment. 
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Transmission lines form a major part of a power system. 
They are exposed to different environmental conditions and 
are spread over a large geographical area. As a_ result, 
lines experience more faults than other system components, 
and so a line protection scheme is essential to any electric 
power network. The types of line faults are mainly 
cia s Sishilied tasssshontaici reuitwand: open (cirewiitine Shortuciirouit 
1S mainly due to insulation failure, human error, lightning, 
wind damage and many others[1]. Open circuit may occur for 
a variety of reasons, including broken conductors and 
maloperation of circuit breakers. The detection, location 
and removal of line faults in the shortest time possible is 


of the utmost importance in the design of line protection. 


The fast and accurate fault location of short-circuited 
transmission lines iS always a prime concern of every power 
utility from economical and safety considerations. A good 
protection scheme will enhance the availability of a power 


system. 


A comprehensive survey of various fault locating 
methods has been given by Stringfield, Marihart and 
Stevens[2]. In this reference, the authors have broadly 
classified the methods as travelling-visual-inspection type 
and measurements-from-fixed-reference-point type. The 
latter type is far superior to the former as less time and 
labor are required to locate a fault. Also, environmental 


conditions may affect the speed in locating the fault for 
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the former type. However, an example of the latter such as 
pulse-reflection type, needs a large amount of expensive 
hardware or well qualified operators to interpret fault 


oscillographic records. 


Recent developments in digital technology have resulted 
in a decrease in cost and an increase in processing’ speed. 
It facilitates the introduction of on-line protection of 
electric power systems on a substation level. The use of 
micro/minicomputers in faut location offers greater 
flexibility and accuracy. In the past few years, new 
short-circuit fault location algorithms based on modern 
techniques have been developed. These methods employ 
measured electrical quantities at one or more points on the 
transmission line. Some of the methods make _ use of 
distributed-parameter line models, and the fault distance is 
determined by analysing the current and voltage travelling 
waves[3][4]. The accuracy of such an approach is affected 
by the system parameters and network configurations. The 
lumped-parameter model algorithm is often used in cases 
where the complex impedance between the faulty point and the 
measuring site 1S calculated. The location can thus be 
obtained by assuming that the line impedance is proportional 
to line length. Thi sem vdea fas Gpopullar Tan ther tield) of 
digital distance relaying. However, the drawback of the 
lumped-parameter approach is that fault resistance and shunt 
Capacitance of long transmission lines introduce error in 


the Fault distance measurement. Dynamic parameter 
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estimation[5] has also been employed in algorithms’ for 
locating faults on power transmission lines. Using this 
algorithm, a set of non-linear equation is formed and solved 
through an iterative approach. In this approach, the 
configuration of the Sequence reference model is determined 


by the type of fault present. 


The block diagram shown in Figure 1.1 outlines’ the 
major components associated with an on-line fault locator. 
The current and potential transformers provide low level 
Signals to the low-pass analog filter. These signals 
correspond to the transmission line currents and voltages. 
The low-pass filter prevents high harmonics and noise 
present in the signals from affecting the fault location 
algorithm. Alternatively, the analog low-pass filter can be 
replaced by a digital one which 1S incorporated into the 
fault location algorithm. This decreases the amount of 
hardware required by the fault locator. The analog to 
digital converter produces the corresponding digital 
information. Fault detection and identification routines 
can be incorporated into the fault distance measurement 
routine to increase the capability of the fault locator. 
Phersoutputs of “the ~Wocator j:provides sinformation for the 


protective relaying system and for off-line fault analysis. 


The object of this thesis is to implement and test an 
existing fault location algorithm based on current and 


voltage measurements at the sending end of a power 
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transmission line. No special line coupling is required 
between the two ends of the line. The method employs least 
error squares curve fitting technique proposed by Sachdev 
and Baribeau[6]. The fault distance is determined by 
assuming that the line reactance is proportional to line 


length between the meaSuring and the fault points[7][8]. 


During the last few years, extensive research has’_ been 
undertaken in the development of fault location algorithms 
using meaSurements of electrical Quantities at some 
observation points. Chapter 2 discusses some of the 


algorithms proposed in the past. 


Chapter 3 presents the algorithm suggested in Reference 
6. The algorithm determines the fundamental current and 
voltage phasors of the post-fault transient waveforms. The 
waveforms include fundamental, third harmonic and decaying 
dc components. A complex impedance is calculated from the 
extracted fundamental phasors and fault distance is obtained 


based on the reactive component of the calculated impedance. 


Digital simulations of transmission lines for testing 
the proposed algorithm are presented in Chapter 4. The 
models provide both pre-fault and post-fault information to 
the weraulLteerlocation aigorithm “in an {off-line mode... The 
development of digital filters for suppressing unwanted high 
harmonics is also briefly discussed in this chapter. The 
cutoff frequencies for the filters are mainly determined by 


the equation of condition derived in Chapter 3. 
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Chapter 5 describes the implementation of the proposed 
algorithm using a PDP-11 microcomputer present in the Power 
Research Laboratory of the University of Alberta. The fault 
location algorithm is tested using transmission line models 


and digital filters which are presented in Chapter 4. 


The performance of the proposed fault location 
algorithm is presented in Chapter 6. The time delay and 
accuracy of the algorithm are investigated using differennt 
order filters. The errors contributed by the presence of 
fault resistance and line shunt capacitance are discussed. 
Certain conclusions pertaining to the validity of the 
proposed fault location algorithm are put forward in the 


final chapter. 
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CHAPTER II 


A review of previously developed fault location algorithms 


Ze ledntroduction 


New digital fault location algorithms for power 
transmission lines have been developed in the past. Most of 
the algorithms are based on meaSured electrical quantities 
at some observation point on the line during the first few 
cycles of the fault occurrence. This approach can _ be 
subdivided into the following four major categories: 

(1) Algorithms based on wave or telegraph equations. 

(ii) Newton-Raphson based algorithms. 

(iii) Reference model based algorithm. 

(iv) Algorithms employing the calculation of complex 


impedance. 


This chapter reviews some of the previously published 
fault distance measurement methods, which belong to one of 
the above four categories. The last category is often used 
in transmission line distance protection since the fault 
impedance loci can be determined based on electrical 


measurements at the relay location. 


2.2 Algorithms based on wave or telegraph equations 


The use of travelling wave phenomena for fault location 
has been presented in some previously published literature. 


With this approach, the fault location is deduced from _ the 
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post-fault current and voltage travelling waves. The waves 


are modelled either by the wave or the telegraph equations. 


Kohlas[3] presented an accurate fault location 
algorithm for power transmission lines. In his approach, 
high frequency transients, due to the occurrence of short 
Circuit Pauley, are incorporated into the distributed 
parameter mathematical model. The meaSured data from one 
end of the line enables the instantaneous profile of the 
line voltage to be obtained using the wave or telegraph 
equation. The fault location is determined from where 
either the voltage profile vanishes or its magnitude reaches 
a minimum. A major drawback of this approach is that the 
effect of a non-zero fault resistance is neglected. In 
addition, the complexity and accuracy of the solution to the 
travelling wave equation largely depends on the overall 


system configuration. 


Another algorithm based on using wave equations, for 
line models, was suggested by Vitins[4]. In this method, 
the fault was located by determining the time delay between 
two quantities at the measurement site. These two measured 
Quantities are associated with the travelling waves. The 
presence of the fault resistance has been taken into account 
in the development of the algorithm. However, a lossless 
transmission line has been assumed and computational effort 


is rather large in analysing the wave phenomenon. 


ar; aigales: 


7 saw moty | Sgr He: et ” Whaat ud ee a, “8 


be pred) Fi 


ade red wiv Cen wail 7a) aodinia hd lisa age han wa 


itd tally ee ston vas, ehis 2a beta 


& 
at 
= 
~ 
| 
— 
2 
ae 
=a. 
a 


vid Pity POs. teen ‘ak ‘ bt hte 128 ta bie: iy 
ot GP Me INK Oe, aie ty yostiotm panne ats ogi 


imeve eit ae Re aga we ‘ope ASL2 mie, svew eaenies 


“ei 
" co 


192 -yenoky aA ween tan dv ie Bvt ao agptR vegtiona 
paseaa wtae.: At Lahiee ay, ud osifes peut baw abso 
chek 2 esc yates din ony erbriaina. te we benea0s ee es 
Baie eeqa owt sige ne! ‘e Spadreess aid a pantry,» 
att done. el teva ate done, busdzvagien shtay peer 
duatos eared abated need aint aaagietzes, phyed me tes 
jaatdeo!, & or ae <a pi ontt a0) om ' 


s1othe Leno seauqnes bia 


2.3 Newton-Raphson based algorithm 


A fault location problem can be formulated into a set 
of nonlinear equations which is solved using the 
Newton-Raphson technique. The Newton-Raphson method is a 
well known iterative method for solving nonlinear systems of 
equations and it has been successfully applied to many power 


system problems[9]. 


Westlin and Bubenko[10] proposed an iterative approach 
to fault location for power transmission lines based on 
Sampled voltage and current at the sending end. A set of 
nonlinear equations is formulated by using a_ suitable 
Thevenin equivalent model for the receiving end of the line. 
Each type of fault results ina different set of equations. 
The fault distance, receiving end current and fault 
resistance can be determined from the solution of equations 
using the Newton-Raphson method. A low-pass filter is 
required to eliminate any high frequency component present 
in the post-fault transient waveforms. The algorithm is 
readily applied to situations where the fault resistance 
cannot be neglected. The disadvantage of this approach is 
that the Thevenin equivalent circuit at the receiving end of 
the line must be known prior to the calculation of the fault 
distance. Also, the iterative process is time consuming and 
hence extensive computational time is needed to locate the 


fault accurately. 
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Anenews types on stavult “location algorithm. susing the 
Fourier transformation method was developed by Takagi, 
Yamakoshi, Baba, Uemura and Sakaguchi[11]. ene eth S 
approach, nonlinear algebraic equations, which contain the 
unknown variable corresponding to the fault distance, are 
formulated by applying the law of superposition to the 
faulted transmission line. Fourier analysis is used to 
extract the fundamental frequency components’ from the 
measurements available at the sending end of the line. The 
extracted components are required by the Newton-Raphson 
iterative program in solving the set of nonlinear equations. 
No communication channels are required between the two ends 
of the line, and the algorithm is able to locate the fault 
accurately without being affected by non-zero fault 
resistance. The disadvantage of this method is that an 
unrealistic lossless line is assumed. In addition, the type 
of source impedance at the sending and receiving end must be 


carefully studied to ensure the algorithm is applicable. 


2.4 Reference model based algorithm 


Recently, the fault location problem was treated as a 
dynamic system parameter eStimation problem by Richards and 
Tan Size |: In their approach, a lumped-parameter 
mathematical model was introduced to compare its response 
with that of the real physical system during the occurrence 
Cie ae reurt. The fault location, fault resistance and 


transformer saturation parameters were computed by varying 
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the model's parameters until an adequate match is obtained 
between the real system and the model response. The 
presence of a fault resistance does not affect the accuracy 
of the algorithm. Due to the asymmetrical nature of a 
fault, a Symmetrical component transformation[9] was used to 
formulate the reference mathematical model. Different types 
of fault result in different configurations of the reference 
model, which increases the complexity of the analysis. 
Also, the sequence reference model can be very complicated 
in the presence of simultaneous faults. In addition, a time 
consuming iterative approach is employed in estimating the 
parameters and hence the proposed algorithm is not suitable 


for on-line applications. 


2.5 Algorithms employing the calculation of complex 
impedance 


A digital distance relay used for line protection 
directly calculates complex impedance between some 
observation point to the faulted point. Ths aeCOnGe Die cor 
impedance calculation is easily applied to the fault 
location problem since the computed impedance provides’ the 


fault distance information. 


Mann and Morrison[13] proposed an algorithm such that 
the peak magnitude of a sinusoid can be determined from its 
value and its rate of change at any arbitrary sampling 


instant. The impedance as seen from the meaSurement site is 
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equal to the ratio of the voltage and current sinusoidals at 
the measuring location. The values of the derivatives at 
any instant are approximated by the difference equations. 
The effects of the presence of a fault resistance and the 
decaying dc components alter the accuracy of the calculated 
impedance. Also, a very small sampling time is required to 


have an accurate approximation for the difference equations. 


A slightly different approach was suggested by 
Gilcrest, Rockefeller and Urden[14]. In this approach, the 
authors used the first and second derivates to calculate the 
peak magnitude of the sinusoidal waveforms instead of using 
the instantaneous values and the first derivatives as shown 
in Reference 13. Three sets of sampled current and voltage 
values are required by the algorithm to determine th complex 
impedance of the faulted line. The benefit of this approach 
is that the decaying dc components present in the samples do 
not affect the accuracy of the algorithm. However, the high 
frequency components introduce errors into the final result 


and thus the fault location is subject to inaccuracies. 


The truncated Fourier Series approach was also used by 
Phadka, Hlibka and Ibrahim[15] in digital relaying problems. 
The algorithm requires samples of current and voltage 
Signals over one-half cycle of the fundamental frequency. 
The samples are then correlated with reference sine and 
cosine functions to extract the fundamental components in 


rectangular form. From these extracted voltage and current 
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components, the phasor impedance can be determined and hence 
the fault can be located. The response for this algorithm 
is fast due to its short data window, nevertheless its 
accuracy 1S not as good as could be since it is highly 
susceptible to high frequency components and decaying dc 


offsets. 


The behaviour of a lumped-parameter transmission line 
model in which the shunt capacitance is neglected can be 
described by the following differential equation at any 


instant of time. 
Vir@rete+ 1(da/at) (2.0 


The above equation is valid for both pre-fault and 
post-fault conditions depending on the choice and 
combination of sampled voltages and currents (v and i). The 
numerical solution of equation (2.1) corresponds to _ the 
series resistance (r) and series inductance (1) of the given 
line. Breingan, Gallen and Chen[16] proposed an algorithm in 
which equation (2.1) is solved at three instances of time to 
calculate the two unknowns r and 1. The advantage of this 
approach is that the decaying dc component is implicitly 
taken into account by the series R-L model of the 
transmission line. The validity of this method is only 
justified for short transmission lines in which the_- shunt 


Capacitance can be neglected. 
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2.6 Conclusion 


Previously published literature related to fault 
location problems has been briefly discussed in this 
chapter. The algorithm based on complex impedance 
calculation is readily adapted to the field of digital line 
protection. This approach is straightforward and requires 
the least amount of computational effort as compared to 
other algorithms. It 1S Suitable for practical situations 
where the fast location of faults are desired. The major 
drawback of this approach is that its accuracy suffers from 
the presence of either fault resistance, decaying dc 
components, shunt capacitance or any combination of these 


factors. 
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CHAPTER III 


Fault Location Algorithm 


3.1 Introduction 


Algorithms using different approaches in solving fault 
location problems have been briefly presented in the last 
chapter. The suitability and limitation of each algorithm 
have also been concisely discussed. The choice of a 
Suitable algorithm that compromises between accuracy and 


speed is a difficult task. 


In this Ghaprer #ijameshort-cireuityetault = locat zon 
algorithm is presented. The method is based on the least 
error squares curves fitting technique proposed in Reference 
6: This curve fitting technique, together with its 
relationship to the pseudoinverse of a matrix, is also 


discussed. 


3.2 Least error squares curve fitting technique and the 
pseudoinverse of a matrix 


Assuming a series of values "a" together with their 
corresponding values "r" are given, it 1s often required to 
find a functional relationship between the two variables r 
and a. Therefore, a set of unknown parameters, which 
describes the relationship between r and a, _ has to be 
calculated. The relationship can be of any function of a, 


inciudingyta straroht. “line, “as ~parabola ‘or acircle. An 
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equation of condition, which describes the known nature of 
tne relatvon or theycurve. to, which sthes,.data’ahasi to be 
fitted, is formulated. Usually, the equation of condition is 


written in the form shown in equation (3.1)[6]: 


Pape Ant xe tae x ft 4. toa ep ey ot an Xn (acai) 


Where oh =e] 22) ba on 
Qvipva+st> Qni- Cm_are the known variables 


X1, eee, X, Are the unknown parameters 


Equation | (Senljameconsisces = Of ~n “unknowns x7 to x, .. In 
general, at least n values of r_ together with the 
corresponding values of a are required to solve for the 
unknown parameters. If n values of r are given, n equations 


of the type (3.1) are obtained and they are shown as below: 
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The above n equations can be written in matrix form, thus: 


ry; aii a12 cece Atn-1 Qin X4 
T2 a21 a22 see Gln 3 QA2n X2 
e = e e e (3:.3) 
Fn 4 Oni en 32 22sec One tines (an - 3 ni i ni- 3 
ln Ant an2 eee Ann-1 Ann Xn 
4 / 
ha R = AX (3.4) 


where, Ris an x 1 column vector 
A is an x n square matrix 


Xe Sia RakepecColuumMaVvector 


5 ho . . . / ° 
If matrix A 1S non-Singular, i1tS inverse A™' exists and 


the following equation is obtained: 
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2 X (A)-'R (365)) 


From equation (3.5), the unknown vector X is solved. Hence, 
the unknown parameters of the equation of condition, which 
describes the relationship between variables r and a are 


completely determined. 


In practice, the number of known variables r is usually 
greater than the number of unknowns =. Under such 
circumstances, the number of equations in (3.2) is larger 


than n and the system of equations is said to be 
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overdetermined[17]. Equation (3.3) is modified as: 


ry a4 ai2 eee | Gd Hires Ain Xj 
T2 a21 a22 ese GQ2n-1 G2n X2 
Um-1 Q@m-11 G@m-12 eee Am-in-1 G@m-in Xn-1 
Um ami am2 eee Qmn-1 Qmn Xn 
or R = AX (3.6) 


where, m>on 
Ris am x 1 column vector 
A is a rectangular matrix of order m xn 


51S aynuxet Column vector 


Equation (3.6) corresponds to an inconsistent system of 
linear algebraic equations. The inverse of A is not defined 
Since A is ye Square matrix. Normally, there is no exact 
solution for the vector X such that AX = R[17]. The least 
error squares technique has to be employed in solving the 


given curve fitting problem. 


The following error-Square expression is obtained for 


the it* equation of matrix equation (3.6): 
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The sum of all the error squares becomes, 
Mm 
Len = e,? + e,? sy eee eB ok a 3 ak 


Hence it follows that, 
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where e;? is the error- square associated with the it 
equation of (3.6). Differentiating equation (3.8) with 


respect to x; yields the following equation, 
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Equating equation (3.9) to zero yields, 
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Rearranging equation (3.10) gives the following normal 


equation[6], 


m m 
dr iais = 2ai1?ky + ai 18 2X2 tee ee tote 
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If equation (3.8) is differentiated with respect to x2 and 
the resulting equation equated to zero and the various terms 


rearranged, a second normal equation is formulated: 
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The above procedures are repeated for x3; to x, and thus 


th 


n normal equations are derived. The n normal equation is 


shown in equation (3.13), 
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whe least . erron squares approach, yields. n normal 
equations from the original m equations. From equations 
Koei) (3. 12)ganedets. 13), ~ thee wormal equations can. be 


expressed in matrix form as: 


Lr jai 1 Lai.’ La j 18; 2 aie sO aMatalin = teod.s aan X 1 
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Or R = AX (3.14) 
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1S an x n Square matrix 
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avltthe summations are from i = 1 toem 


The least error squares process can be summarized in 
the following steps: 

(i) Obtain the error squared expression for each r. 
(ii) Sum all the error square terms. 
(iii) Differentiate the sum of error squares with 
respect to each unknown x. 
(iv) Equate each result from step (iii) to zero and 
rearrange the resulting equation. 


(v) Obtain the final matrix equation R = AX. 
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The definition of matrices A and R, and the product of 
the transpose of A and R gives the following result, 
aii a21 eoe Am-11 ami ry 
a12 a22 eee Am-12 Am2 Y2 
AR = e e e e e 
QAin-1 QA2n-1 eee Qm-i1n-1 Omn-1 Um-1 
Qin QA2n eee Am-in Amn Um 
This can be written in compact form, 
Zrjaii 
Zrjai2 
A'R = 4 
PE eel Wr ieee | 
Zriain 
or: ATR =R (3.15) 
where the summations are taken from i = 1 to m. 
Substituting equation (3.15) into (3.14) yields, 
A’R = AX (3.16) 


(nxm)(Ms1) (nan)(nx4) 


Substituting equation (3.6) into (3.16) gives, 
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OLE Pane = A 32317) 


The least error squares curve fitting problem is now readily 
solved by substituting equations ene .15), and.(.3.407). into 


equation (3.14). The result is, 
A'R = A'AX (3.18) 


TheswnknOWNe. Vector Kes ifng.equation (3.8 ).. is...solved by 
premultiplying both sides of equation (3.18) with the 
inverse of A'‘A, provided that the inverse exists. The final 


equation is given by, 


(puariy ath Ree AA LALA x 


(CX aA. R (2510) 


Ory, x 


(meh ita 1s athemlet tm pseudoinverse of matrix Alo]. Matrix ,A 
must) “have wa, srank? ‘equales tomen. to ensure that A’ A, -is 


non-singular; A'A is invertible[18]. 


In summary, the least error squares curve fitting 
solution is obtained by first finding the pseudoinverse of 
the rectangular matrix A, then the matrix R is premultiplied 
by the pseudoinverse to determine the unknown vector X. AS a 
result, the least error squares problem is resolved into a 
matrix algebraic problem involving the calculation of the 


pseudoinverse. 
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3.3 Development of the equation of condition for the _ fault 
location algorithm 


The least error Squares curve fitting technique and its 
relationship with the pseudoinverse of a matrix have been 
presented in Section 3.2. This technique was suggested in 
Reference 6 as a tool for digital relaying of transmission 
lines. The same approach can also be applied to fault 
location problems. In order to apply the curve fitting 
approach, equations of condition for both the post-fault 


voltage and current waveforms have to be formulated. 


In the presence of a short-circuit fault, current and 
voltage waveforms at the sending end of a transmission line 


are represented by equations (3.20) and (3.21) respectively: 


sek es © + LG, sintkot.+ af) (3772.0) 


v = K*’e-t/7 + FGjsin(kwt + af) 39219) 


where, i is the instantaneous current at time t 
v is the instantaneous voltage at time t 
7 is the time constant of the given system 
K' and KY are the magnitudes of the dc offset of the 
current and voltage signals respectively 
Gi and Gy are the magnitudes of the kt* harmonic of 
the current and voltage signals respectively 
ai and af are the phase angles of the kt harmonic of 
the current and voltage signals respectively 


w is the fundamental frequency of the given system 
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It is assumed that all the fifth and higher harmonics 
present in the waveforms have been eliminated by either a 
low-pass analog or digital filter. Also, the post-fault 
current and voltage signals usually do not consist of even 
harmonic components. AS a result, equation (3.20) is 
Simplified and the resulting equation is shown as (3.22) by 


expanding the exponential term (e-*’7) as a power series. 


eee oat) perenne TO Nem ty / ag leat) 


+ Gisin(wt + ai) + Gisin(3wt + a}) C3222) 


where ! denotes factorial 


The two sinusoidal terms present in equation (3.22) can be 


expanded, thus yielding the following equations: 


Gusimiot + ai) = GiSannot) COS (arp) 


+ Gicos(wt)sin(a}) (3-723) 
and Gisin(3wt + aj) = Gisin(3wt)cos(a3) 
+ Gicos(3wt)sin(a3) (324) 


Substituting equations (3.23) and (3.24) into equation 
(3.22) and approximating the power series of equation (3.22) 
by only the first three terms, equation (3.22) is modified 


into: 


Pos Ki Kit re + K tt 727? 
+ Gisin(wt)cos(a!) + Gicos(wt)sin(a};) 


+ Gisinicut) cos(a,).+ Gicos(3wt)sin(e) (3625) 
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The above equation can be transferred into an equation of 
condition as shown in Reference 6 by using the substitutions 
given in equation (3.27). The resultant equation of 


Condition for the current waveform is Shown below: 


Vera) (o)emetes Ct) x Fay erat a, (t) x. 
+oag(t)xeot ae(t)x¢ + az(tux? C5 25) 

and een ey al (t)ee m4 

X32 = Gicos(a}) az(t) = sin(wt) 

ge NG Keane lt) ali t eacosk wt) 

x4 = G3cos(a}) ai(t) = sin(3wt) KSPIZI79) 

x Mie=2'G YsaniCals ) antGtym= cos Gat) 

xd = -Ki/r eiGre=) tt 

xi = K!/2r? ae CE )ia= ke 


The values of x's are unknown and a's are known since w and 


t are predetermined. 


Proceeding in the same manner, an equation of condition 


for the voltage waveform is obtained and is written as: 


vy Petar’ a acct xa Hralvtt eer thant) xy 


tea SoCt ) omer eth) Get ea Gh) Kis, (3.28) 
where, x) 2) Ke anit) = 1 
XS = Gicosila:) as.(te) Ses in Cote) 
xy = GYsin(a}) a(t) =scos(ot) 
X\= Glcos lay) a, (tee usan (sot) (3.229) 


xalenGs sain las) a¥(t) = cos(3wt) 
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x¥ = -K’/r Salty = 
Kye Shay ote a7(t) = t? 
From.equations (3.27) and (3.29) ; 
ay(t) = aj(t) (3.30) 


Where, ) = 1,2, .s.cnmee 


For simplicity's sake, the superscripts associated with all 


the a's are left out in the subsequent analysis. 


The analog current and voltage signals are converted 
intossdigital \foum pasesthes input, “to. thes fault.location 
algorithm. Each sample of data results in an equation of 
condition either in the form of equation (3.26) or (3.28). 
The number of unknowns, x, of each equation is seven, which 
implies that at least seven samples are required to compute 
the unknowns. In general, a larger number of samples can _ be 
used but each additional sample increases the window size 


and more time is required before the fault can be located. 


At mtheesd’ «Sampling, sinstant) y~uequations, (3.26) and 


(3.28) are rewritten as the following two equations: 


NAS A Ey re Se mes epee 
+ agsX$ + AuggXG + Ag7X7 (3733.1) 
Vv 
and Vor lagikuetmadek cet agsXine toga kn 
Vv 
+ AgexceteagexXs £1ae7X7 (3733:2)) 


where the subscript d denotes the d*tt sampling instant. 
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If m samples of current and voltage are acquired, m 
equations of (3.31) and (3.32) result and these are written 


in matrix form as? 
AX’ « I (3233) 
and AXGL=hy (3.34) 


where, A iS am x 7 matrix 
xX ona s satrem7. x 1 column vectors 


I and V are m x 1 column vectors. 


I and V are known column vectors Since they are consist of 
input current and voltage samples. X' and X” are the unknown 
vectors which are solved from equations (3.33) and (3.34). A 
is specified since the fundamental frequency and _ sampling 


rate are predetermined. 


Usually, the number of samples, m, iS greater’ than 
seven. Under such circumstances, the matrix A becomes 
rectangular in nature. Equations (3.33) and (3.34) are 
solved by applying the concept of least error squares curve 
fitting which was briefly discussed in Section 3.2. The 
unknown values of x are obtained by premultiplying the 


equations with the pseudoinverse of A. 


In the proposed algorithm, the matrix A is determined 
by the sampling rate, the time instant when t is equal to 
zero, and the fundamental frequency of the given system. In 


general, the fundamental frequency is fixed and the 
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remaining two parameters have to be arbitrarily chosen to 


obtain the desired matrix A. 


The first parameter to be considered is the sampling 
rate. The highest frequencey present in the equation of 
condition determines the lower limit of the sampling rate. 
Theoretically, the sampling frequency must be at least twice 
the largest frequency present in the signal to avoid 
aliasing. From equations (3.26) and (3.28), the sampling 
rate has to be at least 360 Hz since the equations consist 
of up to third harmonic components (360 Hz). In practice, a 
Sampling rate larger than 360 Hz is desired to ensure that 
the aliasing effect is totally eliminated from the digitized 
values. If the sampling rate is too high, the determinant of 
matrix A’A becomes small and hence the inverse becomes 
large. A fast sampling rate is undesirable since the 
elements of the matrix (A‘'A)-' have to be multiplied with 
the sampled values and any noise present in the signal would 
be amplified. The sampling rate is determined by considering 


the above factors. 


Another parameter which has to be studied is the time 
reference[6] (the instant when t is equal to zero). 
Different time references result in different matrices A and 
these can affect the implementation of the algorithm. The 
authors of Reference 6 tried three different time references 
such that the time is set at zero at the instant of either 


the first, second or third sample. 
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In order to select the best equation of condition for 
relaying applications, different combinations of sampling 
rate, time reference and sampling window size have been 
Studied in Reference 6. After considering the elements of 
the pseudoinverse of A, the authors decided to select the 
following parameters which result in the pseudoinverse 


Matrix shown in Table 3.1: 


Sampling rate = 720 Hz. 
Total number of samples = 9. 
Time reference = time set to zero at the instant of the 


second sample. 


The above combination of parameters has the following 
advantages: 

(i) The sampling rate is chosen at 720 Hz which is a 
multiple of the fundamental frequency. This eases the 
computational implementation of the algorithm. 

(ii) Nine samples at 720 Hz sampling rate are required 
by the algorithm. This corresponds to three-quarters of 
a cycle at the fundamental frequency which is well 
within the operating time of a protection system (the 
time between the instant of fault inception and the 
time when the circuit breaker trips). 

(iii) Considering the pseudoinverse matrix shown in 
Table 3.1, the values of the elements of the second and 
third rows are symmetrical with respect to the centre 


Sample. Also, one element in the third row is repeated 
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four times and another element is zero. 
These all contribute to the reduction of computation in the 


proposed fault location algorithm. 


3.4 Calculation of complex impedance and determination of 


fault distance 


As mentioned earlier, the unknown vectors X' and X” of 
equations (3.33) and (3.34) are solved by premultiplying 
both equations with the pseudoinverse of A. The resultant 


equations become: 


Reale ex (3735) 
(Tam) (mrt) (711) 

and pCO NBS (35 oy) 
(7am)(mxt) (741) 


where A*[17] is the pseudoinverse of A and is defined as 
(Aan) Ae 


m is equal to nine. 


The second and third elements of the unknown vectors X' and 
XY consist of information regarding the components of 
fundamental current and voltage phasors. By multiplying the 
elements of the second row of A* with the nine samples of 
phase current, the component x2 is obtained. Similarly, x3 
is solved by determining the products of the third row's 
elements of A* with the nine sampled phase currents. The two 
unknowns x} and x3} are solved in the same manner aS x3 and 
x} but the sampled values now correspond to phase to neutral 


voltage. As a result, the following four components are 
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obtained: 


fe ly COS (a) Gers 
X= 1, sintea,) (3338) 
¥76=) ve.cos (ay (3.39) 
X3 = vp Sin(a!) (3.40) 


where, ip 1s the peak value of the fundamental current 
component 
a; is the phase angle of the fundamental current 
component 
Vp is the peak value of the fundamental voltage 
component 
a¥ is the phase angle of the fundamental voltage 


component 


The resistive component R, of the line apparent complex 
impedance as seen from the sending end of the line at any 


instant t can be expressed as: 


R, 


Re {¥] : PRR a eee), 


ys ve [cos(a¥)cos(aji) + sin(a})sin(aj)] 


1p 


ee vp cos(ay)ipcos(a;) + vesin(ay)ipsin(a};) 
Leas ee) 
1p 


p =-Vecos(at)ipcos(a;) + ve Sin(ay)ipsin(aj) (3.41) 
L tp COS a) tp Sl Dato 
Substituting equations (3.37) to (3.40) into equation (3.41) 


yields the final expression of R, [6], 
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Applying a similar procedure, the reactive component X of 
the apparent impedance is shown to be: 


x, = %9Xz — X2x) (3.43) 


Xoxs fox oe 


The calculated R. and xX. from equations (3.42) and 
(3.43) correspond to the resistive and reactive components 
of the apparent impedance of the transmission line 
respectively. In the presence of zero-resistance 
short-circuit fault, the calculated R. is exactly equal to 
the line resistance seen from the fault location algorithm 
to the faulted point. However, fault resistance is 
incorporated into the final R. if the fault resistance 


cannot be neglected. 


Assuming line length is proportional to the line 
reactance, the calculated reactive component of the apparent 
impedance using equation (3.43) is used to locate the 
short-circuit fault. Fault distance ys is determined by the 
following relationship: 


Po aey 4 (3.44) 
Ys = x, 


where, ys is the distance from the sending end of the line 
to the faulted point 
y: is the total length of the line 


Xs is the apparent line reactance from the sending 
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end to the faulted point 


X, is the total reactance of the line. 


y: and X, are known values, y; can be determined once X; is 
calculated. The fault is located using the reactance ratio 
as shown in equation (3.44), and the presence of ground 
resistance does not affect the final fault distance 
measurement. However, the accuracy of the algorithm is 
affected by line shunt capacitance and this effect will be 


Studied in Chapter 6. 


3.5 Special features in using least error squares _ curve 
Fitting approach in fault location 


In the last two sections, the development of the 
proposed fault location algorithm using least error squares 
curve fitting technique has been briefly discussed. Some of 
the advantages of the proposed method are summarized below: 

(i) The freedom in choosing the equation of condition 

which results in the presence of decaying dc component 

and any desired harmonic component. 

(ii) The flexibility in selecting the data window size. 

(iii) Some allowance in choosing the sampling rate. 

(iv) The pseudoinverse can be calculated in an off-line 

mode since the elements of matrix A are all known prior 

to the calculation of the fault distance. 

(v) No iteration is required and hence results in less 


of a computational burden. 
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(vi) The effect of ground resistance to the accuracy of 
the fault location is minimized due to the fact that 
the fault distance measurement is based on the 
reactance ratio. 

(vii) No prior knowledge of source impedance at either 


the sending or receiving end is required. 


3.6 Conclusion 


In this chapter, the least error squares curve fitting 
technique has been briefly studied. It has been shown that 
the pseudoinverse of a matrix provides the solution of the 


least error squares problem. 


The fault location algorithm employing the least error 
squares curve fitting approach has been derived. This 
algorithm is readily implemented by digital computers and it 
has shown that it is well suited to on-line applications. 
The next chapter discusses the various mathematical models 


which are used in the testing of the developed algorithm. 
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CHAPTER IV 


Testing the algorithm using simulated transmission lines 


4.1 Introduction 


A fault location algorithm has been presented in the 
previous chapter. Two simple power systems under both 
pre-fault and post-fault conditions were chosen for study. 
The systems, which were set up in an off-line mode, provide 
transient fault data to the proposed fault location 
algorithm and thus allow the performance of the algorithm to 


be evaluated under different system configurations. 


As mentioned earlier, a low-pass filter is required by 
the algorithm to attenuate high frequency harmonics present 
in the post-fault waveforms. The design of the filter 


affects the performance of the algorithm to a large extent. 


Detection and identification functions can be’ included 
digitally in the algorithm to increase the capability of the 
algorithm. These additional facilities remove the need for a 
communication link between the protective system and the 
fault locating system, thus increasing the portability of 


the fault location algorithm. 


4.2 Modelling of test system 


In testing the performance of the proposed fault 


location approach, two simple power sytems are used. The 
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systems are» Selected 7 in “order ero. represent practical 
Situations. The first one is a Single-phase system 
consisting of a short, aluminium-cable-steel- 
reinforced(ACSR)[19] overhead transmission line. The line is 
mathematically represented by a series equivalent R-L model. 
The second system is a three-phase network consisting of a 
medium length, copper-conductor type overhead transmission 
line. The transmission line is modelled by a_ three-pi 


equivalent circuit.[20] 


4.2.1 Single-phase power system 


A single-line diagram of the single-phase system is 
shown in Figure 4.1. The system, operating under normal 
condition, consists of a generator G1 feeding into a load L1 
through transformers and transmission lines. 111 1s the line 
monitored by the fault location algorithm. The data of the 


components is shown in Appendix A. 


The components of the system are modelled using 
equivalent impedances. Using 126//3 KV as the base voltage 
and 100 MVA as the base power, the impedances are converted 
to per unit base. The steady-state equivalent circuit of the 
system, in which all the impedances are in per unit values, 
is'given in Figlre (4.22 The transmrssron=-line= lit 1s 
modelled by a lumped series R-L network as shown in Figure 
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Figure 4.1 Single-phase test system (pre-fault). 


J° j0-4318 16-3152 j16-9498 
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Figure 4.2 Equivalent circuit of the single-phase test system (pre- 


fault): 


4 1 


The steady-state pre-fault current and voltage phasors 
of the sending end of the line 1i1 are denoted by 1° and vV° 
respectively. From Figure 4.2, the expression of I° and V° 
are written as: 


ica 1.0 Carls) 
J0.4318 + 16.3152 + 316.9488 + 4.8764 + 53.3808 


VES = 1 Ome 1° Oe 43.18 (4.2) 


The presence of a short-circuit ground fault on the 
line 111 yields the faulted system shown in Figure 4.3. The 
Parameter ye corresponds to the distance measured ~ from the 
Sendang, end of the, “fine ‘lit <togethe faulted point. The 
equivalent circuit of the faulted system is shown in Figure 
4.4, The circuit is used to solve for the transient response 


of both if and v* immediately after the fault has occurred. 


In Figure ih iv and vf correspond to the 


instantaneous values. r{ and I. are the’ Series resistance 
and inductance respectively of the faulted line between the 
sending end and the faulted point. Their values are 
dependant upon the fault location since they are directly 


Droportionalitosthne linewlength. © is the generated emf “and 


1Sso Sinusolaal function. 
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Figure 4.3 Single-phase test system (post-fault). 


is 


Figure 4.4 Equivalent circuit of the single-phase post-fault 
test system. 
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The following time domain differential equations are written 


by considering the equivalent model shown in Figure 4.4. 


dt is he 
i.e. dij eae tain (2°35) 
at y+ i 
. § 
and, wi ge rhies if eee (4.4) 


From equations (4.3) and (4.4), the unknowns if and v* are 


solved numerically for each small time step using Euler's 


integration method. The quantities e, l,, ri and i are all 


known prior to the calculation. Initial conditions for if and 


Vv’ are given by the pre-fault steady-state values of I° and 


V° respectively. 


4.2.2 Three-phase power system 


The second test system is almost identical to the 
system described in Section 4.2.1 except now each element 
consists of three components, the a-, b- and c-phase. Line 
lil is replaced by a three-phase 133.5 mile-long overhead 
copper conductor transmission line. The three conductors are 
spaced 15.8 feet apart. Transmission line 1i2 is a 
three-phase 5 mile-long overhead copper conductor line with 
conductors spaced 4 feet apart. For simplicity, no mutual 


coupling is assumed between the three phases. 
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In the normal operating mode, the system can be assumed 
to be in a balanced three-phase state. Figure 4.5 is the 


one-line diagram of the balanced system. 


The single-phase equivalent circuit of the three-phase 
system is shown in Figure 4.6. The various components are 
modelled by their equivalent impedances. The values given in 
Figure 4.6 are in per unit with 126/73 KV and 100 MVA chosen 


as the base voltage and base power respectively. 


The transmission line 1i1, which is monitored by the 
fault location algorithm, is represented by three "pi" 
sections. Each "pi" section corresponds to a 44.5 mile-long 
transmission line, and its line parameters are given by 


equations (4.5) and (4.6) with operating frequency of 60 Hz. 


0.06 + 30.2355 pu (4.5) 


Zx (series impedance) 


Zx (shunt reactance) -355.7 pu (4.6) 


The single-phase equivalent model is used to analyse the 
balanced three-phase system. The distribution of currents 
and voltages is calculated from the network equations. 
Hence, the pre-fault current (1°) and pre-fault voltage (v°) 


are solved. 


If a three-phase ground fault occurs accidentally on 


the line lil, a single-phase post-fault mathematical model 
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Figure 4.6 Single-phase equivalent circuit of the three-phase 


system (pre-fault). 
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bSeused) £O report the condition. This 1S shown in Figure 
4.7. The fault exists at a location 44.5 miles from the 
sending terminal of the line. The use of the single-phase 
equivalent Circuit to study the post-fault system is 


possible since the fault is a balanced three-phase type. 


Using Kirchhotr se current and voltage laws, the 
following circuit equations can be written from the 
post-fault single-phase equivalent circuit (Figure 4.7) 


immediately after the occurrence of the fault: 


PERS eye by 
i.e. Gite e 2 avs (area) 
Cities 1, 
+ 
and, vt = Til + Se 
re, di: vi = reas (4.8) 
TALL If 
of o¢ 
rile ie cise +h 
ie: vee Obani. (4.9) 
dt cr 


The single-phase equivalent circuit is modified and the 
resultant network is shown in Figure 4.8 if the three-phase 
ground fault involves fault resistance ry. The faulted 
network is simplified by assuming no load current and 
current if flows into ground completely through the fault 


resistance. 
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Figure 4.7 Single-phase equivalent circuit of the test system 


(post-fault). 


Figure 4.8 Single-phase equivalent circuit of the test system with 
fault resistance (post-fault). 
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following equations are obtained: 


f Gan 
at, = Rao} 
€ 1,5 
> ait? Mee—* vi (4.10) 
dt i, 
and, -¢ ae : 
Vata ts + SH + M1, 
2 CaP L ev Ebe Rr lider ret (a. iad 
at 1 
and, A f f 
+4 os of + 12 
” CM oo Be NY (se10) 
che ren 


Comparing equations (4.7) to (4.9) and equations (4.10) 
to (4.12), it is observed that the two sets of equations are 
almost identical except for a slight difference between 
equations (4.8) and (4.11). This similarity between the two 
systems eases the computer implementation in simulating the 


two faulted networks. 


Euler's method is employed to solve the differential 
equations (4.7) to (4.9) or equations (4.10) to (4.12) for 
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the unknowns if and v 
initial conditions for the faulted network is obtained from 


the pre-fault model. 


Proceeding in a similar manner as in the 44.5 


mile-fault case, the transient response immediately after 
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the presence of fault is obtained for three-phase ground 
fault occurring at a location 89 miles or 133.5 miles from 
the sending end of the line. Therefore, the proposed 
algorithm can be tested in conditions under which line shunt 


Capacitance is present and fault location is varied. 


In analysing the faulted system in the presence of an 
asymmetrical fault, the Single-phase-equivalent-circuit 
approach is not applicable. Under such circumstances, all 
the three phases have to be considered in order to derive 
the transient response. This analysis becomes very 
complicated if mutual coupling exists between the three 


phases. 


4.3 Design of low-pass digital filters 


The amplitude characteristic of an ideal low-pass 
digital filter with cutoff frequency fc 1S shown in Figure 


4.9. The characteristic has a "brick wall" shape[21]. 


Methods used in the design of digital filter have been 
developed in the past. The text "Introduction to digital 
filtering"[22] edited by Bogner and Constantinides provides 
a good description of the various design methods. One of the 
methods is the "Direct synthesis of digital filters" which 
is used to design recursive low-pass digital filters. This 
technique is employed in this thesis to develop filters for 
the fault location algorithm. The method has a special 


feature such that the Bilinear Transformation[23], which is 
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Figure 4.9 Amplitude characteristic of an ideal low-pass digital 
filter. 
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required to transform transfer functions from s-domain to 
z-domain, has been implicitly taken into account. The design 
problem is resolved into finding a particular function to 
approximate the ideal "brick wall" characteristic when using 


the direct design method. 


Depending on the degree of approximation, different 
particular functions are required in the design analysis. 
One such function is to Simulate a low-pass digital filter 
corresponding to that of a Butterworth-type low-pass analog 
Pelt ers. This function has a magnitude spectrum which 
decreases monotonically in both passband and stopband. Also, 
the magnitude is decreased by 3 dB at the cutoff frequency 
fc. The design procedures of the direct approach is_ shown 
below[22]: 

(we) The order of the digital filter AN)oas chosen from 

the given design specifications. 

(ii) The pole locations are determined on the z~'-plane 

and are selected those that lie outside the unit circle 

to ensure that a stable filter results. Alternatively, 
the poles can be obtained from the z-plane by choosing 
those that lie within the unit circle. 

(in) VANUN on Orcenezenomex)Stswandalceat location zi = 

a6 

(ay) Ae transtersetunetion in’ termsmeot, the Ze! ~ws0r 
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A low-pass filter has to be included in the proposed 
fault location algorithm as shown is Chapter 3. According to 
the equations of condition given in equations (3.26) and 
(3.28), the following specifications have to be satisfied in 
the design of any desired digital filter: 

(i) Fifth and higher harmonic components present in the 

post-fault transient waveforms are attenuated. 

(ii) The sampling rate of the filter has to be at least 


twice the third harmonic frequency (180 Hz). 


Based.) On the direct approach as described above, three 
different ordered digital filters of the Butterworth-type 
have been designed. The details of each filter is described 
in the following. Appendix B gives the design procedures for 


each derived filter. 
(a) 4th order low-pass filter: 


Cutoff frequency = 203 Hz 

Transition frequency = 285 Hz 

Sampling frequency = 720 Hz 

Attenuation at transition frequency better than 30 dB 


Transfer function (G(zn')) is, 


Glz) 6n07 Bole zeae 
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(b) 8th order low-pass filter: 


Cutoff frequency = 200 Hz 
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Transition frequency = 295 Hz 

Sampling frequency = 1440 Hz 

Attenuation at transition frequency better than 30 GB 

Size on) 05023 (ras eee 
UZ 57 oie ze 91235 1) (eee 2362 eae 
a ra a a tee on Le a re ae a Aa 
CZs te More Oto duis 4/5090 eS az a ees S) 


(4.14) 
(c) 10th order low-pass filter: 


Cutoff frequency = 218 Hz 

Transition frequency = 290 Hz 

Sampling frequency = 1440 Hz 

Attenuation at transition frequency better than 30 dB 
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Zi 15 89282-44802 00)} 

The unity step response of the three low-pass digital 
filters is given in Figure 4.10. It has been observed that a 
time delay is associated with each filter. The delay is 
longer for higher order filters. This increases the time 
required to locate the fault. The three designed filters 
are used in the proposed algorithm and the delay caused by 


each filter will be discussed in Chapter 6. 


GEé 


ny | dpe lig ‘rosie & 
~ T7SHeS,S - 


et 
e 7” ~ 
? y xD f f - ’ 


nana’ 193906 ons ant 


eemlign Pete a art we stay ae 


gh i ’ 
exeslt? aaegeeh ~ ager 


ip 
oa 


sh OPS = decaupeat: tieqw 
=a, <<“ DES » qoneupstt cata shane Zs 
sh 0 vps ‘“aheupes2 ont <4 
ey as ae “a 
| oy aha ‘yhaur et a ett 38 not seunettA 


: ivi prea of we «(318 
1 5 Seyi cog Te S-erttsl = Fae) oS 
i 
ve » | 
rg erester’ © # **etex £ =" Ss)” 
; 
Ta ae ; & 
it é § s8S?2.4 - g Ri 
_ } i? 
ig og ¥oL. eatie gi? *0 92 at océess geste tine 
5 ; ' ; oe , : e - 7 


video nod paal, #3 (tam Sheets ai aevip 


ites att aoe kPt apeg Hatw: fernisoses al barry! 


‘ v 
\ 
i | : + t 
7 t 
ak 


imi >. vel uf ott ane widivogia b : oomeney® fd 7% 12 bled 


ont wae aor pwd chee sexo hit 16bI6 sedgid § 


Hol (ge k ey age ai. sdhves sis ‘etal: o2 bs as las JP 
res ye} 5 amar ayer 


nt ae 


ro 


54 


“B494TTS Te{totp ssed-moT 9a14} 04} jo asuodsoa 


ANIL 


yB° Eg ° c2 * Ta° B 


Japio u4t 


Japsto ig 


8 e 
T 
(BT) 4™9T KS ee 
(8) 4>31F3 
Cv) AFIT FS 
qndus deze yl 


YSLI4 IWALISTO G3INSIS30 
JHL 340 
JSNOdS3Y d3LlS =— 29 °T 


¥TUN OT *H eansty 


=2<eOZe-DOwW 


i ae a 


ae ae Pie 
. ese 


Si) 


4.4 Fault detection 


In the implementation of the fault location algorithm, 
two options are available at the start of the fault distance 
calculation. The first is to design a fault detection 
algorithm incorporated with the locating algorithm. The 
detection algorithm initiates the faulted point 
determination routine once the occurrence of a fault is 
detected. The second option is to attempt to locate a fault 
continuously. The latter option has the advantage that no 
detection algorithm is required. However, the disadvantage 
is that meaningless fault locations result from data 
collected from pre-fault state, or from data acquired partly 
between pre-fault and post-fault states. in; adaiztaon, 
unnecessary calculations are needed during the pre-fault 
State since no distinction is made between the pre-fault and 
post-fault states. From the above considerations, a 
detection routine is developed and applied to the fault 


location algorithm. 


The detection algorithm is based on comparing latest 
samples of electrical quantities of all the phases with the 
corresponding samples in the previous cycle. The algorithm 
is operated once a complete cycle of data has been received. 
The data can be either phase current or phase to neutral 
voltage or both. In order to increase the sensitivity of 
the routine, both current and voltage are used in the 


comparison. Also, three consecutive samples are employed in 
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each comparison to accommodate the possibility of spikes 
present in the current and voltage waveforms. Fault 
inception is only acknowledged if all the latest three 
consecutive samples differ from their corresponding 
one-cycle-ago samples by a value which exceeds' the 
prespecified threshold tolerance. The tolerance is an input 
to the detection routine and can be altered by the choice of 


the user. 


4.5 Fault identification 


Logathes proposedwtault) glocatbion « algorithms. a alot. wot 
multiplications and divisions are involved in sethe 
calculation of complex impedance. For a three-phase system, 
Six impedances are calculated to locate all types of 
Shopt-ciurciui tay faults, The S1x impedances are the 
phase-to-phase and phase-to-ground impedances. As a result, 
a lot of computational time is_ required. This greatly 


increases the time needed to locate a fault. 


In order to increase the operating speed of the 
algorithm, it is necessary to reduce the number of impedance 
calculations. High speed is essential for on-line 
applications. One of the solutions is to include a fault 
identification algorithm. This additional routine is 
initiated by the detection algorithm and is used to 
determine the type of fault. Once the fault is identified, 


a suitable choice of current and voltage combination is 
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selected to calculate the appropriate impedance. From this 
calculated impedance, the faulted point is located. 
Moreover, the identification routine can provide a report 


regarding the type of fault for off-line analysis. 


Fault identification algorithms were proposed in 
previously published articles. Reference 16 presents a 


method which belongs to this category. 


4.6 Choice of data for locating fault 


Once the fault type has been identified, a suitable 
combination of current and voltage is required to determine 
the impedance of the faulted phase or phases. For a 
Single-phase-to-ground fault, a compensation technique is 
required to eliminate the sound-phase effect[24] in the 
calculation of impedance. In the case of phase fault, 
"delta"[16] current and voltage are used. Fault resistance 
is included in the impedance calculation if fault resistance 
cannot be neglected. This effect affects the accuracy of 


fault locating. 


If phase "a" of a well-transposed line is subject to a 
single-phase-to-ground fault, equation (4.16)[24] is used to 
calculate the complex impedance. The equation has to be 


modified for the untransposed case[24]. 
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where, Van is the phase to neutral voltage of phase a at 
the sending end of the line 
Ian iS the phase current in phase a at the sending 
end of the line 
I,n 1S the phase current in phase b at the sending 
end of the line 
Ien 1S the phase current in phase c at the _ sending 
end of the line 
K is equal to |2Z,/2Z,| 
Z, 1s the self-impedance of the line 
Zm iS the mutual-impedance of the line 
ZaL is the faulty phase apparent self-impedance as 


seen from the sending end to the faulted point. 


From equation (4.16), the faulted phase to neutral voltage 
and faulted phase current together with a correct proportion 
of non-faulted phase currents are used as the input signals 
for the fault location algorithm. The compensation scheme 
provides a higher accuracy in the final result. Equations 
similar to that of (Cas Sy, are derived for a 


Single-phase-to-ground fault on b- and c~-phase. 


The impedance of a phase-a-to-b fault, phase-a-to-b 
ground fault or three-phase fault is calculated by equation 
(4,17). Similar equations are used for phase faults 


occurring between b and c phases or between a and c phases. 
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where, Van, is the phase to neutral voltage of phase a at the 
sending end of the line 
V_n 18 the phase to neutral voltage of phase b at the 
sending end of the line 
Ign is the phase current in phase a at the sending 
end of the line 
I,n 18 the phase current in phase b at the _ sending 
end of the line 
Z+. 1S the apparent positive sequence line impedance 


between the sending end and the faulted point. 


In considering equation (4.17), the appropriate choices of 
current difference and voltage difference are used for 


faults relating to two or three phases. 


4.7 Discussion 


Two simple power system models and three 
Butterworth-type low-pass digital filters, which are used in 
the testing of the proposed fault location algorithm, are 
presented in this chapter. It has been shown that the 
inclusion of fault detection and fault identification 
routines can increase the operating speed and the capability 
of the algorithm. Compensation by non-faulted phase 


currents reduces the effect due to mutual coupling and hence 
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the accuracy of the algorithm is improved. 


In the implementation of the proposed algorithm, no 
compensation technique is employed in the testing program 
Since the three-phase test model does not consist of any 
mutual coupling between the phases. Also no. fault 
identification routine is included because the type of fault 
is already known during the stage when the faulted model is 
Simulated. However, if the algorithm is applied to 
practical field data, the above two missing routines have to 
be incorporated with the location algorithm to complete the 


algorithm. 
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CHAPTER V 


Computer Implementation 


5.1 Introduction 


In Chapter 3, it has been mentioned that the fault 
location algorithm involving pseudoinverse matrix is readily 
implemented into computer programs. Two programs have been 
developed using "C" programming language[25] based on this 
approach. The simple transmission line models, which were 
presented in Chapter 4, are used to test the proposed fault 
location algorithm. The relationship between the various 
developed programs are outlined in the block diagram as 


shown in Figure 5.1. 


The programs are implemented on a PDP-11 microcomputer 
available in the Power Research Laboratory. The computer 
uses UNIX[26][27] as its operating system and supports high 
level languages such as C and FORTRAN 77. The attractive 
programming environment of the UNIX system and the powerful 


Strength of C language ease the development of programs. 


The developed programs, which are mentioned in Figure 
Sualeameare=) discussed in this chapter. = Thes logic, of each 
programeris > brieflyiapresented: sin ti termShecOtm, general: and 
detailed flow charts. The flow chartS are essential in 


understanding the organisations of the programs. 
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( PSEDO.C ) 


PSEUDOINVERSE 
MATRIX 


PROGRAM POST - PROGRAM_TO 
IMPLEMENT FAULT CALCULATE FAULT LOCATION 
MODEL, FILTER | FILTERED IMPEDANCE AND 
AND AND APPARENT 
MORE FAULT LOCATION 


(MODEL.C) ( LOC.C) 


NOTE : The names in parenthesis correspond to the names of the appropriate 
programs . 


Figure 5.1 Block diagram showing relationship between developed 
programs. 
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59.2 Program "model.c" 


The program "“model.c" is used to eemadace the simple 
transmission line system under both pre-fault and post-fault 
conditions. The test model is either a single-phase or 
three-phase system. In both cases, a single-phase equivalent 
circuit is used to obtain the transient response of the 
faulted system as discussed in the previous chapter. This is 
feasible for the three-phase system since it is subject to 
symmetrical three-phase-to-ground fault. In addition to the 
model simulation, a fault detection routine and a low-pass 


Gigital filter are also implemented in the program. 


5.2.1 Transmission line model routine 


The program starts by declaring the dimensions of the 
variouS arrays and setting the system parameters. The 
Strategy for detecting fault inception is an input to the 
program and is varied according to the user's choice. The 
pre-fault current and voltage phasors, which are obtained 
from the pre-fault steady state model, are required by the 
routine to calculate the instantaneous current and voltage 


prior to the occurrence of a fault. 


The differential equations of the faulted system are 
formulated according to the location of the fault. For a 
non-zero resistance ground fault, the fault resistance has 
to be considered in deriving the differential equations as 


shown in Section 4.2. The equations are solved using Euler's 
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method for the unknown electrical quantities of the sending 
end of the transmission line. The pre-fault and post-fault 
instantaneous current and voltage were printed on a LA120 
Pine printer and piotted ona HP7225 X-vY plotter. Pigure 5,2 
is the general flow chart of the routine which simulates the 


mathematical model of the test system. 


5.2.2 Fault detection routine 


Following the routine of the transmission line system, 
a fault detection routine is implemented into the program. 
The general flow chart of the detection routine is given in 
Figure 5.3. If the current sample lies within the first 
cycle of data, the routine compares the three latest 
Successive samples with a maximum value to decide whether a 
fault has occurred. Otherwise, the routine compares’ the 
three latest consecutive samples with those of one cycle 
ago. In the latter case, the routine concludes a fault has 
occurred if all the three differences exceed a prespecified 


threshold value. 


Immediately after a fault has been detected, one and 
three quarter cycles of pre-fault data and three cycles of 
post-fault data are stored. The pre-fault quantities are 
used to determine the type of fault nya a Pauls 
identification routine is implemented into the program. Only 
three complete cycles of post-fault data are available to 


locate the fault. This simulates a practical Situation in 
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START 
DECLARE DIMENSIONS 
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where 1 is the latest current sample 
SET (absolute value). 
CURRENT V is the latest voltage sample 
tines (absolute value). 


lmex is the maximum allowable in- 
Stantaneous current. 
Vmex is the maximum allowable in- 


SET stantaneous voltage. 
COUNT = 0 tp is the one-cycle-ago current 
Sample(absclute value). 


Vp is the one-cycle-ago voltage 
sample(absolute value). 


ae mek Itele is the current threshold. 
eee eer Veele is the voltage threshold. 


INCREMENT 
CURRENT SAMFLE 


SET 
COUN. = 
COUNT +1 


CALCJLATE THE 

ABSOLUTE VALUE OF 
EACH 

LATEST SANPLE 


. IS . 
[i - tp] > taote 
lv = Vel » Vee 


Figure 5.3 General flow chart of fault detection routine (con'd@). 
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FRINT INFORMATION 
REGARDING 
FAULT 


IS 
FAULT 
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WITHIN 


FIRST 1¢ 
CYCLE 


STORE 2¢ CYCLE OF 
PREFAULT DATA AND 
3 CYCLES OF POST- 
FAULT DATA 


STORE ALL PREFAUL 
DATA AND 3 CYCLES 
OF POSTFAULT DATA 


SET 1=0 AND V=1.0 
FOR DATA RECEIVING 3 
CYCLES AFTER FAULT 

INCEPTION 


PRINT TIME 
RESPONSE 


Pigure 5.3 General flow chart of fault detection routine. 
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which the circuit breaker of the line begins to operate 
three cycles after the instant of fault inception. Once the 
breaker is tripped, the phase current vanishes and this is 


what the routine has simulated. 


5.2.3 Digital filter routine 


The low-pass digital filter routine, which is shown in 
Figure 5.4, is to suppress predetermined high harmonic 
components present in the digitally simulated line data. 
Difference equations of the three developed filters are 
formulated from their corresponding transfer functions. 
These transfer functions have been shown in Chapter 4. The 
difference equations are used to obtain the filtered 
response. The filtered post-fault values are stored ina 
data file which is read by the fault location program to 
determine the faulted point. If it is desired, the filtered 
response can be plotted on a plotter and is compared with 
the unfiltered response. From Figure 5.4, it has been 
noticed that a change of sampling frequency is required as 
the sampling rate is different between the transmission line 


Simulation and the filter simulation. 


5.3 Program "psedo.c" 


In Chapter 3, it has. been mentioned that the 
pseudoinverse of a matrix provides the solution of the fault 


location problem. The pseudoinverse A*, which is present in 
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END 


Figure 5.4 General flow chart of low-pass filter routine. 
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equations (3.35) and (3.36), is calculated from the program 
"psedo.c". The resultant pseudoinverse matrix has already 
been shown in Table 3.1. The matrix is used by the program 
elocec” to calculate: the post-fault “current and voltage 
phasors. The general flow chart of the program "psedo.c" is 


given in Figure 5.5. 


The program starts by formulating the matrix A. The 
elements of A are selected according to the desired 
combination of parameters given in Section 3.3. Matrix 
operators such as transpose, multiplication and inversion 
are all used in the program to calculate the pseudoinverse 
Matrix ({A'’A)~‘'A*)d -The-final—solution is output to a data 
file which is read by the program "loc.c". The advantage in 
uSing two different programs to implement the pseudoinverse 
calculating routine and the fault locating routine is that 
the former routine is only executed once and there is no 
need to execute it every time when the latter routine is 


called. 


5.4 Program "loc.c" 


The impedance calculating and fault distance measuring 
routines are implemented into the program "loc.c". The 
program is to »iocate™ the fault of “the “faulted model 


Simulated by the program "model.c". 


Immediately after the programs "model.c" and "psedo.c" 


have been executed, the program "loc.c" is ready to be 
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Figure 5-5 General flow chart of the program to caiculate 
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called. Figure 5.6 is the general flow chart of the program 
"loc.c". The inputs to the program are the pseudoinverse 
matrix, order of the employed filter, total line length, 
entire line impedance, and total number of post-fault 
Samples along with the respective current and voltage 
values. Equations (3.42) and (3.43) are used to calculate 
the apparent impedance seen between the sending end of the 
line to the faulted point. The location of fault is then 
determined from the computed reactance ratio using equation 
(3.44). The final result is plotted on a X-Y plotter and if 
desired, printed on a line printer. From Figure 5.6, it has 
been observed that the program continuously locates the 
fault until the last set of nine samples has been employed. 
If a 4‘ ordersfiiter iseused in the program “model:c"™, nine 
Successive samples are required LOT the impedance 
Calculation. In the: case ofa6' = or 10°” ordermetylter, every 
second sample is skipped to collect nine consecutive 
samples. This is necesSary to simulate the sampling 
frequency of the fault location algorithm since its sampling 


frequency is half that of the digital filter. 


5.5 Conclusion 


This chapter has outlined the programs which have been 
developed in order to implement and test the proposed fault 
location algorithm. Flow charts are given to explain the 
details of the appropriate routines. The programs are 


programmed on a PDP microcomputer to study the 


“Asuna ‘sale “faz01 


tiusi~ 330g 16 een ine 
psilov hae Baers out 


me (ae. 8 ane 
siz to Bee gothaye ady + naa eA. nese san 
tte? @1 ‘fy 4 96 , sods oad “it \gntog’ besiue a 

iyoupe Oiiny oftey ba ea slashed deel ota 0 2 Be nt kare 


— 


nite wHTIO%9 F-H we, ne ati io le as ‘fluee? tenid sar fi i? a 


statuotes of bese 4x6 008 ‘7 


| — ey 
1wph4, Gos .. “tesmi sa emit a no berniig 


(are bes 
+ oe: tt ted a 
evolan . Be p02 sesh: one ran us ht 


ita ,%. Labor" maorpady eda pawimt 2 $2684. teaae ee 
| a ons — i =e 
e,relegus qj7 «= 307 Tie bogey Sand |  aslghse evies 9920 
etizt rchira Or ao” “OS. 3e eens aia a .eobtstun: 34. 
‘ y _ 1, ; 

tuspengs sale shell’ oF, wea! te 22 sigue 
peliquesa fj. stelamie of eit ai a2 


eri ce need aos el ro nitiyouls ios tee com shied oid , 


rettie Lavieds of¢ de ast iLed ate 


need evan doitw rawsoen 98 senttwe eit 29: igi | eidt 
airs a ligud ea" s 
a a 


7m 
ee : 


“en ate: Pr 


ei 


START 


READ IN 
POSTFAULT FILTERED 
SAMPLES 


READ IN A* 
LINE LENGTH 
LINE IMPEDANCE 


READ IN 
ORDER OF 
FILTER n 


SET CURRENT SAMPLE 
TO FIRST SAMPLE 


RETRIEVE NEXT 
8 SAMFLES 
FROM MEMORY 


RETRIEVE Ne es 
EVEN SAMFLES 
FROM MEMCRY 


CALCULATE SUMMATICN 
USING LATEST 
NINE SAMPLES 
AND A* 


Figure 5.6 General flow chart of program "“loc.c" (contd). 
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Figure 5.6 General flow chart of program "loc.c". 
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performance of the algorithm using differently ordered 
filters as applied to differently faulted systems. These 


studies will be presented in the next chapter. 


CHAPTER VI 


Performance of fault location algorithm 


6.1 Introduction 


The proposed fault location algorithm as described in 
the previous chapter was implemented using computer 
programs. Simple power transmission systems, which are 
Ourlined “in Chapter 4, were simulated to obtain data for 
testing one fault location algorithms. The following 
sections in this chapter will present the performance of one 
of the algorithms under different conditions. These 
conditions include employment of different order low-pass 
filters, variation of fault location, presence of fault 
resistance and line shunt reactance in the _ faulted 


transmission line model. 


The time response of the current and voltage signals, 
which are measured at the sending end of the studied power 
line, will also be given inthis chapter. Once a fault 
occurs along the line, transients are produced in both the 
sending-end current and voltage. These transients are 
affected by both the location of fault and the instant of 
fault inception. In testing the algorithm, the fault is 
initiated at an instant such that the generated line to 
neutral voltage of the faulted phase has a phase angle of 
240°. If necessary, the instant of fault inception can be 


easily varied by. changing the appropriate input of the 
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program "model.c" as mentioned in Chapter 5. 


Ideally, once a fault is initiated, both the collected 
current and voltage samples can be used in the calculations 
of impedance and fault distance. However, there is a delay 
of response associated with each filter and this causes a 
delay in acquiring samples by the fault location algorithm 
to perform the appropriate calculations. For example, if the 
third post-fault sample together with its eight successive 
Samples are used as input data instead of using the first 
nine post-fault samples, a delay of two sampling intervals 
results for the corresponding fault distance estimation. 
This delay in locating a fault affects the performance of 
the algorithm to a great extent since a desired fault 
location algorithm has to locate the fault in the shortest 
possible time. In this chapter, the delay in acquiring data 
for calculating the fault location due to different order 
filters is studied. Figures of calculated fault distance 
versus delay are given. All these reSulting distances are 
obtained from post-fault data with delay of at least 
0.0125 s. If samples with a delay of less than 0.0125 s are 
used, the computed fault distances are very much different 
from the true distance and these meaningless results are not 
shown in the appropriate figures. Since the circuit breaker 
of the line is assumed to be tripped three cycles after the 
fault inception, the calculated locations shown in the 
figures of fault location versus delay are all based on 


samples collected within the first three post-fault cycles. 
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Calculated impedances under different testing 
conditions are also presented in this chapter. The 
impedances are plotted on R-X plots. In each plot, the 
impedances are calculated from different sets of samples by 
varying the delay in acquiring samples until the first three 
cycles of post-fault data have all been used. Similar to the 
estimated locations shown in the appropriate figures, the 
calculated impedances given in the R-X planes are all based 


on samples with not less than 0.0125 s of delay. 


6.2 Single-phase test model 


The single-phase power transmission system, which is 
outlined in Section’ ©4.2.1,) 4S Used ing wstudying the 
performance of the proposed algorithm. The algorithm is 
applied to a 30 mile-long overhead transmission line. A 
short-circuit ground fault is applied at a distance 10 miles 
from the sending end of the line. No resistance is involved 
in the short-circuit path. Theoretically, the impedance 
between the sending end and the faulted point is 5.44 + 
j 5.65 pu. Since there is no fault resistance and no line 
shunt reactance present in the faulted model, a high degree 


of accuracy is expected in the fault distance measurement. 


The transient response of the faulted system before and 
after filtering are given in Figures 6.1 through 6.6. These 
figures are obtained by employing different ordered low-pass 


filters to suppress undesired harmonics from reaching the 
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algorithm. From these figures, it has been observed that no 
high frequency component is present in the unfiltered 
post-fault waveform. In this special case, a low-pass filter 
is not necessary in testing the algorithm. However, a filter 
is still employed as an added safeguard. Considering the 
post-fault current response, the peak magnitude is almost 
three times that of the pre-fault value. This large current 
causes overheating to the system components if the fault is 


not cleared in the shortest possible time. 


Figures 6.7 and 6.8 show graphs of calculated fault 
distance versus delay and the corresponding calculated 
impedances respectively when a 4*th order filter is 
incorporated with the fault location algorithm. It has been 
observed that if samples with a delay of at least 0.0275 s 
are used, the corresponding estimated fault location is very 
close to the true value of 10 miles. Longer delay increases 
the accuracy of the final result. The calculated impedances 
converge to the true impedance of 5.44 + j 5.65 pu as_ shown 


in Figure 6.8 if longer delay is permitted. 


The performance of the algorithm, when an 8" order 
filter is employed, is given in Figures 6.9 and 6.10. If a 
delay of at least 0.029 s is allowed in acquiring samples by 
the algorithm, both the calculated impedance and estimated 
fault distance are almost equal to their respective true 


values. 
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mne «solutions yor) the fault Miocation and impedance 
calculations are plotted on Figures Guid and Oral? 
respectively for the case when a 10** order filter is used 
in the testing. The simulated results compare well with the 
true values if the corresponding calculations are based on 
input Samples collected at least O2032 45 after the 
occurrence of the fault. If samples acquired before 0.032 s 
are used, the algorithm is inaccurate and the resultant 
impedances are scattered over the R-X plane as seen in 


Bagureso. 12: 


6.3 Three-phase test model 


In Section 4.2.2, mathematical models of a simple 
three-phase power transmission system under both pre-fault 
and post-fault conditions are developed. These models are 
used in testing the validity of the proposed fault location 
algorithm. A three-phase URE) mile-long overhead 
transmission line is monitored by the algorithm. The line is 


w 


modelled by three sections of an equivalent "pi CIrCuicte. 
Three-phase ground faults occur along the line at distances 
of 44.5 miles, 89 miles and 133.5 miles from the sending end 
of the line. The effect of non-zero fault resistance on the 
accuracy of the algorithm is also studied. This non-zero 


resistance ground fault occurs at a location 44.5 miles from 


the generating terminal of the line. 
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The unfiltered and filtered responses of “the test 
System pprior to gang after ‘the cecurrencemor amraultware 
given in Figures 6ei@ethrough 6.18. The various Figures are 
obtained by varying the location of fault. The responses 
shown in these Figures correspond to the case when an 8th 
order filter is employed for filtering purposes. Similar 
resultsrare obtaineduusingragd: naorder filter and a 10°* 
Order hilterm=™ whichtMaresegot » Given im this thesis. From 
Figures 6.13 through 6.18, it can be observed that both a 
decaying component and high frequency components are present 
in the pre-filtered post-fault transients. The high 
frequency components are especially noticeable in the 
voltage waveform when the line is subjected to a fault 
occurring at the sending end. The filtered responses, as 
shown in Figures 6.13 through 6.18, illustrate that the high 
frequency components have been suppressed. These indicate 
that the designed filter functions properly according to the 


design specifications. 


The estimated locations and calculated impedances’ that 
result when using the proposed algorithm, together with 
different order filters, are given in Figures 6.19 through 
6.24. These results are obtained when a fault is initiated 
at a distance 44.5 miles from the sending end of the 
monitored line. The measured distances agree favorably with 
the true fault location (44.5 miles) provided that a certain 
delay in acquiring data is permitted. The true impedance is 


0.061 + j 0.236 pu and compares well with the calculated 
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Figure 6.17 Time response of the pre-filtered and post-filtered current signal using 8th order 
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Figure 6.20 Plot of calculated impedances - one pi faulted model and 4th order filter. 
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Figure 6.22 Plot of calculated impedances - one pi faulted model and 8th order filter. 
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Figure 6.23 Plot of fault location versus delay - one pi faulted model and 10th order filter. 
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impedances as seen in Figures 6.20, 6.22 and 6.24. AS seen 
PhomeEagures: 6.19) eee) and 6.23 09a Minamunemdeda yess 
necesSary in order to have the calculated distance converge 
to a reasonable accuracy. This minimum delay varies between 
different ordered filters and is longer for higher ordered 


filters: 


The results of the fault location and impedance 
calculation shown in Figures 6.25 through 6.30 are for the 
case where a fault occurs at a distance 89 miles from the 
generating terminal. The computed locations and impedances 
are within the desired accuracy when compared to the true 
values. Among the three filters, the 10'* order filter gives 
the highest accuracy in locating a fault although a longer 


delay is needed to have the final result converge. 


Figures 6.31 through 6.36 illustrate the performance of 
the algorithm if a fault is present at the receiving end of 
the simulated line. The performance is acceptable for both 
the 8t* order filter case and 10*t* order filter case as 
shown in Figures 6.33 through 6.36. Both the fault location 
and calculated impedance converge to the desired values. 
However, the final result of the 4t* order filter case is 
not satisfactory. The calculated impedances scatter all over 
the R-X plane and the computed locations do not converge to 
the true value. One possible explanation to this undesired 
effect is that the high frequency components present in the 


post-fault waveforms are not completely eliminated by the 


episvios worm 


Naow 4 ad | ask saul ; 


oongoocmt ba “ne seme Posh oie iT) ealunes 
si3 202 ote DE 4 dogdaita, #8. 4 evel ab morte no? 


ont Mowe aelin ay sonatielly ei eu sat. tivet. é saisiw 98 
| | | vs. 
eganebeqm? fee ena taa0% Sesugess ed? . Lea tmx wh ane 1 : 


a 
aa 


#32 aud? +92. as etme teiw ‘gouiv>es beviest “at? ninztw | . 
watit jubye: wge ott .87 tt wits eds prema. aoutew 


Oe “s ee i : 
mecot’ a d¢voielé Sit eens del (ve Wwls8 asta 


agiaenes: 1 bier» (99 sal wradiat heSisan at gale 


opigetigiies s07 sgbiieut ving wo 6 ‘a eomupit 
‘> boo @Aivicdey sta ge Socedty et stent, * ti oistroete # 
Joa 3O2 eldaercoowe BE eanensoi FQ ah vomit bevetunta 
2% an80 sad¢il? ibaa 01 dita eae festit 39670, a. it | 


ms! wan ie 
Hos IeIOL tive? aff sha an. 4 oD eka zetupit: mee nvorde 


| bas 
reufay bealesh ed (ag epreoneh waeruibagad noes eo 
ei eaeD rails dye ae edt 30 hue tenit ote on ak 
I9ve Ite 19273 ana . salle ber oleokaw — ee r 
ot eeeawncs ton v6 anetsmaot besupnen aa i 


‘bei hesbne aids om) noi 


ot od smeneng pibeees p< 


eda ed betsninite” 


oe le 
> 


OMmrm~eaeEIi ZOownwaernor 


mMmozranrmysa 


25 
LOCATION VERSUS DELAY 


~4TH ORDER FILTER 
185 


KEY 
+ loc. of fault 


188 


95 


85 
81 .@15 . 82 . B25 . 83 . 835 . 84 

TIME-S 
Figure 6.25 Plot of fault location versus delay ~ two pi faulted model and 4th order filter. 


PARAMETER IDENTIFICATION (4) 
BY CHANGING 


-So Tt THE INSTANT OF ACQUIRING DATA 


-53 

+ 

+ 

ao | 

+ 

+ 
. 49 + 
iden. impedance 
- 47 
.- 885 . 895 oat e115 wiico 


RESISTANCE 
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Figure 6.27 Plot of fault location versus delay - two pi faulted model and 8th order filter. 
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Figure 6.28 Plot of calculated impedances - two pi faulted model and 8th order filter. 
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Figure 6.30 Plot of calculated impedances - two pi faulted model and 10th order filter. 
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Figure 6.32 Plot of calculated impedances - three pi faulted model and 4th order filter. 
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Figure 6.33 Plot of fault location versus delay - three pi faulted model and 8th order filter. 
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Figure 6.34 Plot of calculated impedances - three pi faulted model and 8th order filter. 
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Figure 6.35 Plot of fault location versus delay - three pi faulted model and 10th order filter. 
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filter, and these unwanted components cause an aliasing 
effect on the digitized signals since the sampling rate of 


Ene digital filter is only 720 Hz. 


In order to study the error introduced by the presence 
of fault resistance, the algorithm is applied to a line 
Subject to non-zero resistance ground fault! occurring ‘at a 
location 44.5 miles from the sending terminal. The 
resistances used are 0.063, 0.126 and 0.189 pu respectively. 
Figure 6.37 gives the plots of location versus delay for the 
three different values of fault resistance. The plots are 
all based on employing an 8t* order filter to suppress 
undesired harmonics. The three curves converge to three 
different values which are quite close to the true fault 
location (44.5 miles) as shown in Figure 6.37. Accuracy in 
locating a fault is higher if the fault involves a smaller 
resistance. This effect is illustrated by the three graphs 


CEigurne 6.37). 


6.4 Conclusion and General Discussion 


The performance of the proposed fault location 
algorithm has been studied in this chapter. From these 
tests, it can be concluded that the algorithm is able to 
locate faults with very good accuracy even in the presence 
of fault resistance and line shunt capacitance. The 
low-pass filter, which is required by the algorithm to 


eliminate high undesirable harmonics, causes a time delay in 
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locating fault. This delay has to be minimized in order to 
increase the availability of power transfer. A low order 
Pricer 1S) suitable ine Shorteningaadelay. timess. but. ithe 
accuracy in fault location is not as high compared to 
Situations where a higher order filter is used. Asa 
result, the choice of filter is a compromise between speed 
and accuracy. Since all tests are based on data obtained 
from digitally simulated simple models, it is always 
desirable to test the algorithm utilizing off-line field 


data. 


One possible future project is to evaluate the 
performance of the proposed algorithm for on-line 
applications. The microcomputer system, which is available 
in the Power Research Laboratory, is well suited for this 
purpose. The system consists of an LSI 11/2 microprocessor 
Whthwesldiotes of wsefuljsiburit—im«x components pidesigned:scfor 
on-line and real-time applications. This complete system is 
known as SNORKEL[28]. The programs, which are outlined in 
Chapter 5, can be easily modified in order to be executed in 
the SNORKEL environment for on-line testing. The modified 
programs together with the SNORKEL system can result in a 


useful on-line fault locating tool. 
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CHAPTER VII 


Summary and Conclusion 


A tast.and:) accurate Jlocations,of faults ton. cower 
transmission lines is essential to minimize the down-time of 
an electrical power system. The availability and security 
of the system is increased and maintenance costs associated 
with locating faults by manual methods are reduced. A 
short-circuit fault location can be determined either by 
manual, visual inspection Or by making electrical 
measurements from some monitoring points on the lines. It 
has been shown that in terms of speed, the latter method is 
more favorable than the former. This is especially true if 
the lines are long and run through inaccessible areas. The 
latter approach is gaining acceptance due to the advantages 
offered by recent progress in both the analog and digital 
technologies. The cost and flexibility of microprocessor- 
based monitors are factors influencing the acceptance of 


automated fault location. 


The method of using electrical measurements in locating 
faults involves either the modelling of transmission lines 
using wave equations or by lumped-parameter representation. 
The wave theory describes a complex phenomenon and is 
difficult to analyse even with the aid of digital computer. 
Among the various fault location approaches based on 
lumped-parameter representation, it has been shown that’ the 


one which employs the calculation of complex impedance is 
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the least complex in terms of computational burden. No 
iterations are required by this method and thus it is 
well-suited for on-line applications. However, one of the 
common drawbacks of this algorithm is that it ignores the 
fault resistance, decaying dc components present in the 
transients, line shunt capacitance or any combination of 
these factors. These factors all contribute to the loss of 


accuracy in the final fault distance measurements. 


An impedance calculating algorithm, proposed by Sachdev 
and Baribeau, iS a promising one which is readily 
implemented into a computer program. In this algorithm, the 
Current and voltage phasors are obtained using curve fitting 
technique. From these phasors, the impedance seen in the 
direction of fault from the sending terminal of the faulted 
line 18S calculated. Hence, the location of the fault is 
determined by assuming that the line reactance is 
proportional to line length. One of the special features of 
this algorithm is that the decaying dc component has_ been 
taken into consideration during the development of the 
algorithm. The effect of fault resistance on the fault 
distance determination iS minimized due to the fact that the 


fault is located on the basis of reactance ratio. 


Mathematical considerations have shown that the fault 
locating problem can be formulated into a rectangular matrix 
form by suitable choices of data window and sampling rate. 


The problem is then solved by a least error squares curve 
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fitting technique in which the pseudoinverse of the 


resulting matrix provides the solution. 


Mathematical models of two simple test systems working 
under both normal and faulty conditions have been derived in 
this project. These models are essential in supplying both 
Prestault "and: post-rault "data “to teste therfaurt location 
avcorirthm. in addition oO unis,  LOW-passe digital. fi1iters, 
which suppress undesirable high frequency components and 
noise, have been developed. These filters are necessary for 


the evaluation of the performance of the algorithm. 


Computer programs were developed to implement’ the 
proposed fault location algorithm in an off-line mode. It 
has been demonstrated in this thesis that the proposed 
algorithm satisfactorily locates faults for the simulated 
transmission lines. The performance of the algorithm is 
acceptable even in the presence of fault resistance and line 
shunt capacitance. However, the algorithm loses its 
accuracy if the post-fault load current cannot be neglected. 
Further studies on the proposed fault location algorithm 
using both off-line and on-line field data are necessary 
with particular attention placed on testing the algorithm 


subject to different fault types. 
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Appendix A 
Power System Parameters 


Data for the system given in Figure 4.1 is outlined in 
this appendix. This data is given in References 19 and 20, 


and is given here for ease of reference. 


(i) Generator G1: 
105 MW capacity at 0.9 pf rated at 22 KV 


18.1 percent reactance on 105/0.9 MVA 


ae short-circuiting Limiting» ceactore 


me Ole 23780 


ovyeecoanstormer Ti: 
Transformation ratio 142/22 
Load carrying capacity of 60 MVA 
8.2 percent reactance on 60 MVA base at 142 KV _ tap 


setting 


(iv) Transmission line lil: 
30 mile-long overhead aluminium-cable-steel-reinforced 
conductor 
Series resistance of 0.5438 Q2/mile 


Series reactance of 0.5650 Q/mile 


(vy) Transformer T2: 
Transformation ratio 126/33 
Load carrying capacity of 20 MVA 


11 percent reactance on 20 MVA base at 126 KV tap 
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setting 


(vi) Transmission line 112: 
5 mile-long overhead copper conductor line 
Total series resistance of 1.5 Q 


Total series reactance of 3.3 2 


(via), LoadaLid * 
Real power of 15 MW 
Reactive power of 8 MVA 


Termanal’ voltage of asi. SaKY 


Appendix B 
Synthesizing Digital Filters 


Three low-pass recursive digital filters were designed 
and employed in the testing of the proposed fault location 
algorithm. The design specifications for the filters are 
given in Section 4.3. The design of tthe recursive filters 
are based on the "direct" method outlined in Reference 22. 
The synthesis technique for each filter that was used in 


this project is briefly outlined in this appendix. 
B.1 4th order filter 


Mines cutoft (ic) transityvon (£1) “andisampling (fs) 
frequencies of the filter are all given in the design 


specifications. Hence, the sampling interval is given by: 
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2.945905 


and, tan(2s*) = tan( 2555) 
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At the transition frequency, the following equation is 
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obtained by considering the amplitude characteristic: 


Mofo] Wate 30 
tan(292) 
hes per £62.999565 
] (2s) 
1.223939 
i.e. Nee oes eee oe 


where n is the order of the required filter. 


Using the calculated value of tan(wcT/2), it can be shown 
elatatrne = poles; which “l1e@ron -the sz =plane, haver ithe 


following x and y coordinates: 


ae eee ee) 
-0.104637 #0.0968 17 
-0.144995 02656432 
ble, Siaiispehets) +1.448533 
a2, LU Sowo ae chcslis) || 7474.2) 


Tt) is »noted “that four “poles ive insidesthe unit circle=of 
the z~'-plane and the other four lie outside the unit 
Circle. Therefore, the four poles which yield a stable 


function are: 


-0.318989 + j31.448533 


and -2.105975 + 33.961229 


The corresponding polynomial due to the stable poles is 


given by: 
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Deve ce + OCI EB OCeae 44058) (zo 0316905 -a 71. e4 eos) 
KZ) +82. 10590" ta S ao 6 ac eect? 10598 (= 4359612 5)) 


SZ +o 0638070) st 2,20) tae tee ez i 2024 +20 51265) 
Thus, the required transfer function is: 
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Gli nodes sia Gee} 


where K iS a gain constant and is given by: 


Be 


1 = K(37§37908)(25.33842) 
6.07802 
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B.2 8th order filter 


The sampling interval T is given by: 


ari 
y= fic 
Paes: 
= Ta40 ° 
a! tan(“3t) = tan(Tea0) 
rT 5014 
and, We T = 2007 
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n 1s obtained from, 


on = —2-999565 
Is (9-75014) 
J\ 0.46631 
erg ae n= 7.26 = 8 


Using the calculated) value of tan(weT/2),.1t) ‘can. be “shown 
that the x and y coordinates of the poles which lie on the 


z '-plane are: 


X y 
Ons6/027 +0208 5336 
Um 20S +0.259994 
0.450890 +0.446794 
ORO 22a +0.653643 
Or LOO 2:7, +0.883342 
1.119040 SieUGO7S 
1.770487 te a2254 
2.584865 +0.600985 


The eight poles which yield a stable function are: 


0.755727 + j0.883342 
1.119040 + j1.108873 
1.770487 + j1.172254 


2.584865 + j0.600985 


Therefore, the polynomial in terms of z~' which corresponds 


to the stable poles is: 
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cee) 7 Pollo zee + 1, Sk 2) Cee. 238420 41254819) 
Rize mae 34 Ozmeuet 4.5080) (zm e b> o1697c> Sat 7.0477) 


Hence, the required transfer function is given by: 
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where K is calculated from: 
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B.3 10th order filter 


From the design specifications, the sampling interval 


is given by: 


fs 
2 1 
1440 ° 
and, Grins 200%) 
tan (23 ) PEO Era 
= O27 53:23 
and, t Ors = o187 
es: ) S5n rae 
Zi 1503 
Hence, n is calculated from: 
2.999565 
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Using the calculated value of tan(weT/2), bs and y 


coordinates of the poles are shown to be: 


wee Zeya 
0.321884 +0.070592 
0.336568 +0.214214 
0.368548 +0.365347 
0.423998 +0.529629 
0.515107 +0.713253 
0.665457 +0.921438 
Geet 73 +1.150664, 
1.368509 +1.356622 
2.114579 +1.345853 
2.964173 +0.650062 


The ten poles which yield a stable transfer function are: 


0.665457 + j0.921432 
0.921173 + 31.150664 
1.368509 + 51.356622 
2.114579 + 51.345853 


2.004 out) J 0.650062 


Therefore, the corresponding polynomial due to the stable 
poles is given by; 
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Hence, the required transfer function is: 
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where K is given by: 
Be 


= Krongeosn nGimsaon a1 somes. NesuO5sou) (4 2e0e2) 


aloe K = 0.032246 
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